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Following an earlier analysis which examined the effect of the self-gravity of pressure on big-bang 
nucleosynthesis (BBN), we explore the effect of pressure's self-gravity on the structure of neutron 
stars. We construct an ad hoc modification of the Tolman-Oppenheimer-Volkoff equation wherein 
pressure's self-gravity is parameterized by a constant, Xi with < x ^ 1- The full general relativistic 
contribution to the gravity of pressure is recovered with x = b an( i i s eliminated when x = 0. This 
formulation is not proposed as an alternative theory of gravity, but is merely used to quantify the 
extent to which the self-gravity of pressure contributes to the structure of dense objects. As can 
be surmised qualitatively, neutron star masses can be quite sensitive to x, with higher values of 
neutron-star mass (by ~20-25%) allowed for smaller values of X- However, for a given equation of 
state, neither the range of neutron star radii nor the radii at fixed central density depend sensitively 
on x- Over the neutron star mass range measured so far, the presence or absence of pressure's 
self-gravity yields a nearly immeasurable change in radius — much smaller than the variations in 
radius due to the uncertainty in the equation of state. In contrast to the result for BBN, we thus 
find that neutron stars are not likely to be useful testbeds for examining the self-gravity of pressure. 
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I. INTRODUCTION 

One of the more fascinating predictions of general rela- 
tivity is that pressure is self gravitating. This is a strictly 
non-Newtonian result, whereby the pressure of a field 
contributes to gravity, increasing the effective density by 
an amount 3P/c 2 for a homogeneous perfect fluid. 

There arc four prominent physical situations in which 
the self-gravity of pressure could potentially lead to a 
measurable effect (three of which involve the expansion 
history of the Universe): (i) the acceleration of the Uni- 
verse during inflation; (ii) the expansion of the Universe 
during the radiation-dominated epoch; (hi) the current 
acceleration of the Universe due to dark energy; and (iv) 
the mass-radius relation for neutron stars. 

In an earlier paper, Rappaport et al. [l[ addressed (ii) 
by setting a constraint on the self-gravity of pressure dur- 
ing the radiation-dominated epoch of big bang nucleosyn- 
thesis (BBN). They introduced an ad hoc multiplicative 
parameter x to the 3P/c 2 contribution in the Einstein 
field equations. For x = 1, the full general relativistic 
contribution to the self-gravity of pressure is retained; 
for x = 0, there is no such contribution. This yielded a 
modified set of Fricdmann-like equations. The contribu- 
tion to (a/a) 2 for a species of matter with equation of 
state (EOS) w = P/pc 2 appears as 

\aj \ l + 3w J o^ 3 ™ V ' 

The gravitational effect of radiation (w — 1/3) is there- 
fore scaled by (1 + x)/2- The absence of presure's self- 
gravity means that the Hubble constant at the time of 
BBN would be smaller by a factor of \/2, potentially 
changing the abundance of light elements. Rappaport et 
al. used a standard BBN code and performed light ele- 
ment abundance calculations for a range of x to quantify 



this effect. When combined with current light element 
observations, the data were shown to be consistent with 
X = 1, and strongly exclude x = 

Our goal here is to see whether similar limits can be set 
by observations of neutron stars. At neutron star densi- 
ties, P ~ (0.1 — 0.6)/oc 2 , suggesting that pressure is high 
enough that a measureable self-gravity effect might ex- 
ist. We introduce a parameterization of self-gravity very 
similar to that used to modify the Friedmann equations. 
The result is a parameterized equation of stellar struc- 
ture, with x = 1 reproducing the Tolman-Oppcnheimcr- 
Volkoff equation, and x — "turning off" the pressure 
self-gravity of that equation. 

Although at fixed central density x nas a significant 
impact on a neutron star's mass, it has very little impact 
on its radius. If measurements were to determine both 
mass and radius, it would still be extremely difficult to 
tell the difference between models with x = and x = 1 
due to the uncertainty in the neutron star EOS. Neu- 
tron stars thus appear to not be very useful for testing 
the self-gravity of pressure. This is not to say that one 
cannot make interesting statements about gravity with 
neutron stars 0; but, in contrast to the situation with 
BBN, the pressure self-gravity aspect cannot be usefully 
tested. The difference between the two cases is simple: 
the EOS of the universe during BBN is well understood, 
but the EOS of neutron stars is not. Testing pressure's 
self gravity is thus degenerate with testing the EOS for 
neutron stars, but is not degenerate during BBN. 



II. STRUCTURE OF NEUTRON STARS 

In general relativity, the structure of spherically sym- 
metric fluid equilibria such as neutron stars is gov- 
erned by the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tion ■ To motivate our heuristic, it is useful to review 
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the derivation of the TOV equation. Further details can 
be found in, for example, Ref. Sec. 23.5. 

We begin by assuming static, spherical symmetry in a 
perfect fluid. The perfect fluid stress-energy tensor is 



= ( fi + P)u, 



(2) 



where p is the fluid's local energy density, P its pressure 
(with c = 1), and is a component of its 4- velocity. In 
equilibrium, the fluid is static so its only non-zero com- 
ponent is ut. The metric follows from the line clement 



ds 2 



3 2 * dt 2 + dr 



2A j„2 , r 2 dQ 2 + r 2 sin 2 ^2 _ (3) 



The condition g^u^u" = — 1 sets u t = — e*. From the 
metric, we compute the Einstein tensor G^ and enforce 
the Einstein field equation (EFE) G^ v = 87rGT pl/ : 



G tt = ^e 2( *- A > (2rf - 1 + e 2A ) = SrrGe 2 * p (4) 
£ (2rf + 1 - e 2A ) 



8nGe 2A P . (5) 



Because of spherical symmetry, the Gee and G<j,$ compo- 
nents do not contain any additional information. 

Local conservation of energy, expressed as V^T' 1 " = 
(where V M denotes the covariant derivative), gives 



dP 

dr 



-(P + P) 



d<I> 

dr 



(6) 



Note that (p + P) plays the role of an inertial mass den- 
sity here. It sets the proportionality for the force a fluid 
element must feel in order not to experience free fall. 

With a little rearranging and a judicious definition, Eq. 
(jl|) describes how the star's gravitational mass accumu- 
lates as a function of radius: 



m(r) = l(l-e- 2A ) 



4tt / p(r')(r'fdr' . (7) 



This result may appear intuitively obvious. Note, how- 
ever, that 4irr 2 dr is not the proper spherical spatial vol- 
ume in the metric ((3]); it lacks a factor of y/g rr - The 
difference between m(r) and a mass m'(r) that includes 
this factor can be interpreted as the gravitational binding 
energy of the star. 

Combining Eqs. ||7J) and (0 into Eq. ^ at last yields 



dP 
dr 



= ~{P + P) 



G[m{r) +3P(47rr 3 /3)] 
r 2 [l - 2Gm(r)/r] 



(8) 



This is the standard form of the TOV equation. Note the 
varying roles P plays here. The p + P term acts as an 
inertial mass density. By contrast, m(r) + 3P(4-7rr 3 /3) 
acts as a gravitational mass. In any metric theory of 
gravity, this mass distinction arises because the dP/dr 
term represents a force accelerating a fluid clement with 
inertial mass density {p + P) away from its geodesic. 

Thus, only the second P term in Eq. (|SJ) represents the 
gravitating component of pressure. We now parameterize 



the contribution of pressure to dP/dr by modifying the 
TOV equation as follows: 



dP 



P) 



G[m{r) + 3Px(47rr 3 /3)] 
r 2 [l - 2Gm(r)/r] 



(9) 



In contrast to the modified Friedmann equations dis- 
cussed in the Introduction, \ — does not correspond 
to a Newtonian limit. The [1 — 2Gm(r)/r) term in the 
denominator is a statement of geometry, and the (p + P) 
term arises from conservation of energy that must be 
valid even in the special relativity limit. 

It is worth emphasizing that Eq. ([9]) is not derived by 
proposing a plausible alternative theory of gravity. Our 
goal is much less ambitious than reformulating gravity; 
we merely want to quantify the extent to which the self- 
gravity of pressure contributes to the structure of dense 
objects. The point of this exercise was to see which of 
the pressure terms in Eq. © acts as a source of gravity, 
and then to appropriately modify that equation in a way 
that tests that term's importance. 



III. MODEL BUILDING, RESULTS, AND 
DISCUSSION 



The first step in building neutron star models is to 
select an EOS, p = p{P). For high densities, we use 
the compilation of Lattimer & Prakash @|, plus two of 
the "mixed" EOSs from Alford et al. [J Our set in- 
cludes a wide range of models, representing several differ- 
ent theoretical approaches, which incorporate both nor- 
mal nuclear matter compositions and more exotic matter 
(quarks, hyperons, etc.). TableUgives an EOS summary, 
listing general properties and giving relevant references. 
The trailing numerical labels are as given in Q . For the 
"ALF" EOSs, the numerical label indicates the transition 
density in multiples of the nuclear saturation density. In 
all cases, we use the Negele-Vautherin (NV) EOS [8[ for 
the lower baryon density range 0.08 > > 6 x 10~ 4 
fm~ 3 , and we use the Baym-Pcthick-Sutherland (BPS) 
EOS @ for n b < 6 x 10~ 4 fm~ 3 . This is similar to the 
low density choices made in Q. 

We construct the neutron star models by numerically 
solving the coupled differential equations 



dp 

dr 
dM 

~~dV 



dp dP 
dP~dr~ 

4:Trr 2 p . 



(10) 
(11) 



The tabulated EOSs were logarithmically interpolated 
and the derivatives evaluated by a simple 3-point La- 
grange interpolation. Integration proceeds until the den- 
sity falls below that of iron, 7.9 g cm -3 . We constructed 
models with both x = and x = 1 f° r each EOS. 

Figure 1 illustrates the effect that varying x has for a 
particular EOS, AP4. We show the mass-radius relation 
for x = 1 (black curve) and x = (red curve). No- 
tice that the radial range spanned by these two models 



3 



Symbol 


Reference 


Approach 


Composition 


AP(3-4) 


Akmal & Pandharipande [10] 


Variational 


np 


MPA 


Miither, Prakash, & Ainsworth fill 


Dirac-Brueckner HF 


np 


MS(l-3) 


Miiller & Serot [12] 


Field Theoretical 


np 


WFF(l-3) 


Wiringa, Fiks, & Fabrocine [13] 


Variational 


np 


ALF(2-3) 


Alford et al. [7] 


Quark Matter 


npQ 


GM(l-3) 


Glendenning & Moszkowski [14] 


Field Theoretical 


npH 


GS(l-2) 


Glendenning & Schaffner-Bielich [15] 


Field Theoretical 


npK 


PCL 


Prakash, Cooke, & Lattimer [16] 


Field Theoretical 


npHQ 



TABLE I: Summary of the EOSs shown in Figure [2] and [3] "Approach" characterizes the theoretical basis. "Composition" is 
characterized by (n - neutrons, p - protons, H - hyperons, K - kaons, Q - quarks). The horizontal line separates "normal" 
and "exotic" equations of state. 



is largely the same; the mass, however, can change sub- 
stantially. The arrows in this figure connect models that 
have equal central densities. This shows that pressure's 
self-gravity can have a ~10 - 30% effect on the mass, but 
barely change the star's radius — at a fixed central den- 
sity. One might have guessed on intuitive grounds that 
an effect of order tens of percent must occur, since for 
these models P/p ~ 0.1 - 0.6. The fact that this change 
is mostly in mass, leaving the radius relatively unaffected, 
was not obvious in advance. Finally, in this regard, we 
note that for a given EOS and at a fixed neutron-star 
mass near the upper mass limit, the radius can vary by 
up to ~15% when \ is reduced to near zero. However, 
the differences in radii at a fixed mass, among the var- 
ious plausible EOSs, are much larger than this, thereby 
making it difficult, at best, to constrain \ 

Figure [5] shows the effect of removing the self-gravity of 
pressure for several illustrative EOSs. The general trend 
seen in Fig. [I] — an increase in mass for a given radius — 
holds for both "normal" (left panel) and "exotic" (right 
panel) compositions. In the most compact stars, \ = 
corresponds to a change of ~ 0.5Af Q over the nominal 
X = 1 model, for a constant radius. 

In the case of BBN, the cosmic abundance of light ele- 
ments puts interesting constraints on x- Might it be pos- 
sible to place similar limits on x from neutron star obser- 
vations? The answer, unfortunately, appears to be "no" . 
Neutron star masses have in some cases been measured 
to accuracies of greater than 1 part in 1000, most no- 
tably PSR 1913+16 [13) j where the pulsar and its com- 
panion have measured masses of 1.441M Q and 1.387Af Q 
respectively. A measurement of the radius of one of these 
objects would place strong constraints on the EOS. Re- 
ferring to Fig. [2l in this mass range and for "normal" 
EOSs, it would take a radius measurement with an accu- 
racy ^300 m in order to, for a given EOS, differentiate 
between the x = an d X = 1 cases. A number of pro- 
posed techniques for measuring neutron star radii have 
been recently summarized in [l8fl . While inventive and 
intriguing, none of these techniques has yet led to robust 
determinations of neutron-star radii. 
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FIG. 1: For a representative EOS (AP4), the mass-radius 
relations for GR (x = 1, black) and without the self-gravity 
of pressure (x = 0, red) are compared. For the x — 1 an d 
X = curves we connect models with equal central density. 
We note that the last two arrows on the left side of the plot 
point from stable x = 1 models to unstable \ — models. 
Thus, there are no stable models where the radius changes 
by more than ~ 5% due to variation in \ f° r this particular 
EOS. 



Figure [3] makes this point even more clearly, including 
a more comprehensive sample (all EOSs listed in Table 
1) and showing the ranges of possible masses and radii as 
shaded regions, bounded by the most extreme EOSs. The 
change in the region between the x = and x = 1 cases 
clearly shows that the absence of gravitating pressure 
acts to stretch possible masses upward, coupled with a 
minimal increase in radius. 

On this plot, the only plausible evidence for a model 
with x 7^ 1 would be a precise measurement of both mass 
and radius that placed a star in a red-shaded region that 
is not also a black-shaded region. For neutron stars with 
M < 2 Mq it would be extremely difficult to distinguish 
between viable EOSs on the one hand, and the effect of 
pressure's self-gravity on the other. For higher masses, 
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ALF3 X = 
GM2 x-1 
GS2 
PCL 




2 Mq < M < 2.7 Mq, a measurement of radius could in 
principle distinguish the effect of self-gravity from EOS 
dependence, at least for those EOSs we have considered 
here. Of course, if any neutron stars are ever found with 
masses exceeding ~2 Mq, the entire set of EOSs would 
certainly be robustly reexamined before considering that 
the equations of general relativity should be modified to 
remove the self-gravity of pressure. 
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FIG. 2: Mass-radius relations for a variety of EOSs. The 
general relativistic result (x = 1) is shown in black, while the 
results in the absence of the self-gravity of pressure (x = 0) are 
shown in red. The left panel shows EOSs of nuclear matter. 
The right panel contains equations with more exotic kinds of 
matter. 
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FIG. 3: Mass-radius relations for a spectrum of EOSs. The 
shaded area fills the region bounded by the most extreme 
M(R) curves. The general relativistic result (x = 1) is shown 
in black, while the results in the absence of the self-gravity of 
pressure (x = 0) are shown in red. 



Given the difficulty of interpreting such a measurement 
as a test of gravity, we conclude that neutron stars (in 
contrast to BBN) are not a good laboratory for examin- 
ing the self-gravity of pressure. 
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